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■ Abstract 

o ; 

, We compare the predictions of the coUinear parton model and the A:7^-factorization 

I approach in J/^' and D* meson photoproduction at HERA energies. It is shown that 

^ ' obtained in the both approaches D* meson spectra over px and rj as well as J/^' meson 

Pt- and z-spectra are very different. The predictions obtained in the fcy-factorization 
approach are agree with the experimental data well. We also predict the nontrivial 
PT-dependence of the the spin parameter a{pT) in the J/^' photoproduction. 
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O . 1 Hard processes in the parton model and /c^- factorization 

g:^ '■ approach 
o 

■ Nowadays, there are two approaches which are used in a study of the charmonia and charmed 

mesons photoproduction at high energies. In the conventional collinear parton model pf] it 
is suggested that hadronic cross section, for example, o"(7p ccX,s), and the relevant 
partonic cross section d'{'yg cc, s) are connected as follows 



O 



{iP ~^ ccX, = J dxG{x, iM^)a{'jg — > cc, s), 



^ . where s = xs, G{x,iJ,'^) is the collinear gluon distribution function in a proton, x is the 

^ I fraction of a proton momentum, /i^ is the typical scale of a hard process. The evolution 

of the gluon distribution G{x,ii^) is described by DGLAP evolution equation [2]. In the 
so-called fc^-factorization approach hadronic and partonic cross sections are related by the 
following condition [SI IH Ej : 

a^^(7P-cX,s) = j^jdklj g$(a:,4,/iV(7^7^-cc,S,4), (2) 

where a{jg* cc, s, k'^) is the cc-pair photoproduction cross section on off mass-shell gluon, 
k'^ = kx^ = — fcy is the gluon virtuality, s = xs — k^, Lp is the azimuthal angle in the 
transverse XOY plane between vector k^ and the fixed OX axis. The unintegrated gluon 
distribution function $(a;, can be related to the conventional gluon distribution by 



xG{x,i/) = j ^{x,k'^,i/)dk'^, (3) 







^Email: saleev@ssu.sainara.ru 
^Email: vasin@ssu.samara.ru 



1 



where $(a:, k^yfi"^) satisfies the BFKL evolution equation [0]. In formulae (2) the four- vector 
of a gluon momentum is presented as follows: 



k = xpN + kr, (4) 

where kx = (0, AJj-jO), pn = {E]\f,0,0, Ej^) is the four-vector of a proton momentum. At 
the a; -C 1 the off mass-shell gluon has dominant longitudinal polarization along the proton 
momentum. Taking into account the gauge invariance of a total amplitude involving virtual 
gluon we can write the polarization four-vector in two different forms: 

e'{k) = (5) 



or 



ef^^k) = (6) 



As it will be shown above formulae (5) and (6) give the equal answers in calculating of 
squared amplitudes under consideration. 

Our calculation in the parton model is down using the GRV LO (?] parameterization 
for the coUinear gluon distribution function G(x, /i^). In the case of the fc^-factorization ap- 
proach we use the following parameterizations for an unintegrated gluon distribution function 
$(a;, k'^.iJ?): JB by Bluemlein |H]; JS by Jung and Salam KMR by Kimber, Martin and 
Ryskin JU]. The detail analysis of the evolution equations lied in a basis of the different 
parameterizations is over our consideration. To compare different parameterizations we have 
plotted their as a function of x at the fixed fc^ and yU^ in Fig.l and as a function of k"^ at 
the fixed x and /x^ in Fig. 2. 

Note, that all parameterizations of an unintegrated gluon distribution function describe 
the data from HERA collider for the structure function F2(x, Q^) well [Hj 1^ ITU]. 

2 meson photoproduction in LO QCD. 

The photoproduction of the D* meson was studding experimentally by HI and ZEUS Col- 
laborations at HERA ep-collider {E^ = 27.5 GeV, En = 820 GeV) Because of the 
large mass of a c-quark usually it is assumed that D* meson production may be described 
in the fragmentation approach ^ ^3]; where 

a(7p -> D*X, P) = J Dc^D* {z, /i^)cx(7p cX, pi = p/z)dz (7) 

and Dc^D*{z, fi^) is the universal fragmentation function of a c-quark into the D* meson 
at the scale fi^ = m'jj + p^. The fraction of the D* produced by a c-quark as measured by 
OPAL Collaboration 



± 

J Dc^D*{z, fx^)dz = 0.222 ± 0.014, 



^The another approach based on recombination scenario was suggested recently in |14| 
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has been used in our LO QCD calculations to normalize the fragmentation function. 
The Peterson pfH] fragmentation function was used as a phenomenological factor: 



c^D 



( 2n _ Ar ^(1 ~ 



In the high energy limit or in the case of a massless quark one has following relation for the 
four- vectors p = zpi, however in the discussed here process the D* meson energy is not so 
large in compare to M^* and the following prescription was used 



p = zpi 



(9) 



together with the mass-shell condition for the c-quark energy and momentum Ef = 01 + m^. 
We have used e = 0.06 as a middle value between two recent fits of D* meson spectra 
in e"'"e~-annihilation, which based on massive charm (e = 0.036) ^H] and massless charm 
(e = 0.116) [T7| calculations. The squared matrix element for the subprocess ■jg* cc after 
summation over a gluon polarization accordingly (6) may be written as follows |H IH} I18j: 



M 



167r^e^asQ; ■ (s + k 



T) 



{t — m1){u — ml) 



2ml 



02 



2, 



u — m^ 



t — m? 



(10) 



where s, i and u are usual Mandelstam variables, 



ai 



mi + pIt mi + p^2T 
, a2 



t 



mi: — u 



Pit and p2T are the transverse momenta of c- and c-quarks, fcy = Pit + P2T- 

Using formulas (5) for a BFKL gluon polarization four-vector we can rewrite (10) in the 
another form: 



IMP 



IGn'^elasa 



[ml — tY{m1 — uy 



m„ 



2ml - 4:mlpf^k^ + m^k^ + Spf^k^ + 



+3k^ + (4m^ + Upi^k^ + 5k^)s - (3m^ - Apf^ - 3k^)s^ + + 



5m" 



imlpfrpk"^ — 2mlk^ — Apfrpk^ — k^ — 12m^s — Af^rpk'^s — 
-kp + 6mls'^ - kls^ - s^\i - Upi^kl - k^ + 3{-2ml + sY ]t' + 



+4:(2ml - sjP - - 4|pit| MA?t| cos(^)(-2m 



- (4 - s - 2i)i{2ml - u) + m^ikl + 3s + Qt) + ml {3k^ + f + 
+kl{As - 2i) - 6si - 6P)) + |pit| cos(2y?) (m^A-^ + kli{2ml - u) - 
-ml{2k^ + s^ + kl{3s + 2t)^ 



:iii 



where (p is the angle between piT and kx- 
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In the last case (11) it is easy to find the parton model limit: 



27r 



where 



and 



\MpM? 



hm /^|M|2= |Mpm|2, 

' 



|fcT|-*0 



^ _ ^ _ {u-ml){i-ml) 2 
Pit ~ P2T ~ ' ""^ci 



(12) 



((m2-t)2(-m2 + s + t)2) 



-2ml + 8ml{s + i) ~ 

-t(s^ + 3sH + AsP + 2P) + m^(s3 + Gs^t + 8P + 4st(3t + u)) 
-m\{lf? + 12P + 4s (4t + u)) 



(13) 



3 meson photoproduction at HERA 

In this part we will compare our results obtained with leading order matrix elements for 
the partonic subprocess — > cc in the conventional parton model as well as in the fcr- 
factorization approach with data from HERA ep-coUider. The data under consideration 
taken by the ZEUS Collaboration ^1]. Inclusive photoproduction of the D*^ mesons has 
been measured for the photon-proton center-of-mass energies in the range 130 <W < 280 
GeV and the photon virtuality <\ GeV^. At low the cross section for ep eD*X 
are related to 7p cross section using the equivalent photon approximation jT^ : 



where f-y/eiu) denotes the photon flux integrated over from the kinematic limit of Q^in = 
mly'^/{l — y) to the upper limit Q^ax = 1 GeV^, y = W'^/s, s = AEnEs, En and Ee are the 
proton and electron energies in the laboratory frame. 
The exact formulas for f^/e{y) is taken from j20j : 



My) = ^ 



a l + il-yf , ^ , 1 



log + 2miy 



The limits of integration over y are y^ax = Wmax/s. In our calculations we used formulas 

min min 

for differential cross section in the following form: 



da{ep eD*X) f r . , . f dz ^ , 2^ i d^p /• -2'^(x, A;|,,/i2) 

- y dyf^i,{y) J —D,^D*{z,fi^) J ^ J dkl x 

^b-.lb-.H IMP 
EiiW^ -2En{Ei~pi,)) WttxW^ ^ ^ 



drjdpj 



The differential cross section as a function of the D* pseudorapidity, which is defined as 
r] = — ln(tg I) , where the polar angle 9 is measured with respect to the proton beam direction, 
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is shown in Fig. 3 where the kinematic ranges for the D* meson transverse momentum are 
2 < pt < 12 GeV, 4 < < 12 GeV, 6 < px < 12 GeV, correspondingly from up to down. 

We see that the results of calculations performed in the collinear parton with LO in 
matrix element need additional i^'-factor [K ^ 2) to describe the data. The value of this 
K-factoT is usual for a heavy quark production cross section in the relevant energy range. 
Opposite the results obtained in fcr-factorization approach are agree with the data well. 
Especially in case of the JB parameterization of a unintegrated gluon distribution function. 
The difference between the predictions obtained with the various parameterizations is about 
factor 2 near the peak of the cross section da/drj. Note that our results coincide with the 
results of the Ref.^] in case of the JB parameterization. 

The pt spectrum of D* meson in photoproduction at |?7| < 1.5 and 130 <W< 280 GeV 
are shown in Fig. 4. The theoretical curve in the fc^-factorization approach obtained using 
the JB parameterization describes data well at the all pt- The curves obtained using KMR 
or JS parameterizations coincide with the data only at the Pt < 8 GeV. The typical value 
of the i^-factor in case of LO parton model calculation is equal 2 at the large pt- 

Our results show that the introducing of a gluon transverse momentum in the frame- 
work of the /cT^-factorization approach increase the D* meson photoproduction cross section 
at the large pt in case of the JS parameterization, as it is predicted for the J/\E' photopro- 
duction [221 122] (see next part of the paper) and at the all relevant values of the px in case 
of KMR or JB parameterizations. 

The dependence of the D* meson production cross section on a total photon-proton 
center-of-mass energy W is shown in Fig. 5. Nowadays the experimental data for da/dW are 
absent. We see that the difference between the results obtained with the various parameter- 
izations of a unintegrated gluon distribution function is about 50%. As well in the case of 
the 77-spectra. At the all pTmm the cross section calculated in the parton model is smaller 
than predictions obtained in the /cr-factorization approach. 

The main uncertainties of our calculation come from the choice of a c-quark mass in the 
partonic matrix elements (10), (11) and (13), and from the choice of a parameter e in the 
Peterson fragmentation function Dc^£,*{z, fi^). 

We see that agreement between the results obtained in the /cT-factorization approach and 
the data on the D* meson photoproduction is better at the large px- This fact shows that 
in the fc^-factorization approach the next to leading order corrections are taken into account 
effectively. 

4 J/ip photoproduction in LO QCD 

It is well known that in the processes of J/ip meson photoproduction on protons at high 
energies the photon-gluon fusion partonic subprocess dominates [21]. In the framework of 
the general factorization approach of QCD the J/ip photoproduction cross section depends 
on the gluon distribution function in a proton, the hard amplitude of cc-pair production as 
well as the mechanism of a creation colorless final state with quantum numbers of the J/ip 
meson. In such a way, we suppose that the soft interactions in the initial state are described 
by introducing a gluon distribution function, the hard partonic amplitude is calculated using 
perturbative theory of QCD at order in ^^(/i^), where ^ ~ nic, and the soft process of the 
cc-pair transition into the J/ip meson is described in nonrelativistic approximation using 
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series in the small parameters as and v (relative velocity of the quarks in the J/ip meson). 
As is said in nonrelativistic QCD (NRQCD) j2H], there are color singlet mechanism, in which 
the cc-pair is hardly produced in the color singlet state, and color octet mechanism, in which 
the cc-pair is produced in the color octet state and at a long distance it transforms into a 
final color singlet state in the soft process. However, as it was shown in papers [23 120] , the 
data from the DESY ep-coUider [27] in the wide region of pt and z may be described well 
in the framework of the color singlet model and the color octet contribution is not needed. 
Based on the above mentioned result we will take into account in our analysis only the color 
singlet model contribution in the J/ip meson photoproduction [25. We consider here the 
role of a proton gluon distribution function in the J/ ip photoproduction in the framework of 
the conventional parton model as well as in the framework of the fc^-factorization approach 

[HllllSI- 

There are six Feynman diagrams (Fig.6) which describe the partonic process 'jg — >■ J/ipg 
at the leading order in as and a. In the framework of the color singlet model and nonrela- 
tivistic approximation the production of the J/ip meson is considered as the production of 
a quark- ant iquark system in the color singlet state with orbital momentum L = and spin 

momentum S = 1. The binding energy and relative momentum of quarks in the J/ip are 

P 



neglected. In such a way M = 2mc and Pc 



Pc 



-, where p is the 4-momentum of the 



J/ip, Pc and Pc are 4-momenta of quark and antiquark. Taking into account the formalism 
of the projection operator |28j the amplitude of the process •yg J/'ipg may be obtained 
from the amplitude of the process jg — > ccg after replacement: 



V%p,)U^{Pc 



(15) 



where i{p) = e^(p)7^, £^j.{p) is a 4- vector of the J/ip polarization, —z= is the color factor, 

\E'(0) is the nonrelativistic meson wave function at the origin. The matrix elements of the 
process 'jg* J/'ipg niay be presented as follows: 

M, = Krt„(gl)e^(g)4(g2)e.(p)Mf^'^^ 



(16) 



Ml 



m: 



Tr 



Tr 



Y{p + M)7^ 



Pc 



{Pc - qi 



qi + mc -p-c 

2^ 



rat 



{Pc 



- q2 + mc p 
- ml 



ut- , nj\ Pc + q2 + mc ^ k-pc + ruc 
7 {p + M)-f''- ■ ^7 ^7^ 



(Pc + q2Y - ml {q- pcY 



mt 



(17) 



a0fiu 



m: 



af3fiu 



-- Tr 
: Tr 
Tr 
: Tr 



7^(p + M)7" 



Pc- qi + mc p 



{Pc - qi] 



mt. 



7 



Pc + mc 

2 ' 



{q -PcY - ml 



k + m, 



{Pc - qf 



c a Pc 

7 



mt 



{q2 + p 



q2 + mc Q 

2r 



mt 



, p Pc + q2 + mc qi- Pc + mc 



{Pc + q2Y - ml (gi - pcf 



mt 



-.1 



Y{p + M)Y 



Pc-k + mc 
{Pc - qY - ml ' (gi - PcY - ml 



p qi- Pc + mc „ 
-.Ttz. -?7 



(19) 
(20) 
(21) 
(22) 
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where qi is the 4-momentum of the photon, q is the 4-momentum of the initial gluon, q2 is 
the 4-momentum of the final gluon, 



K = e^e^^^lW fjab ^ 1 TrfT'^T''], = -, e = Vte, gs = \/4^. 

The summation on the photon, the J/ip meson and final gluon polarizations is carried out 
by covariant formulae: 

^c,{qi)ef^{qi) = -g^p, (23) 

spin 

£a{q2)epiq2) = -gai3, (24) 

spin 

E^.iP)^Ap) = -9,u + ^. (25) 

spin 

In case of the initial BFKL gluon we use the prescription (5). For studing J/ip polarized 
photoproduction we introduce the 4-vector of the longitudinal polarization as follows: 

In the high energy limit of s = 2(gip7v) >> M'^ the polarization 4-vector satisfies usual 
conditions {elBl) = —1, {^lp) = 0. 

Traditionally for a description of charmonium photoproduction processes the invariant 
variable z = (pPn)/ (QiPn) is used. In the rest frame of the proton one has z = E^/E^. In 
the /cT^-factorization approach the differential on and z cross section of the J/ip photo- 
production may be written as follows: 



dp^dz z{l 



JO 27r Jo IQ-nixsY 



The analytical calculation of the |Mp is performed with help of REDUCE package and 
results are saved in the FORTRAN codes as a function of s = {qi + q)'^, t = {p — qi)'^, 
u = {p — qY, p|., and cos{ip). We directly have tested that 

dip 



hm /"-^|M|2 = |Mpm|2, (2J 

:2^0J0 iTT 



tu 



where P5- = — in the IMP and IMpmP is the square of the amplitude in the conventional 

parton model |24j. In the limit of k|. = from formula (27) it is easy to find the differential 
cross section in the parton model, too: 



daP^'\^p ^ J/mX) ^ \MpM?xG{x,p.'') 
dp\dz lQ'K{xsyz{l — z) ' 

However, making calculations in the parton model we use formula (27), where integration 
over and (p is performed numerically, instead of (29). This method fixes the common 
normalization factor for both approaches and gives a direct opportunity to study effects 
connected with virtuality of the initial BFKL gluon in the partonic amplitude. 
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5 J/ijj photoproduction at HERA 



After we fixed the selection of the gluon distribution functions /x^) or <!>( ) there 

are two parameters only, which values determine the common normalization factor of the 
cross section under consideration: \E'(0) and m^. The value of the J /ip meson wave function at 
the origin may be calculated in a potential model or obtained from experimental well known 
decay width T{J/ip — >• /i"*"//"). In our calculation we used the following choice |\E'(0)p = 0.064 
GeV^ which corresponds to NRQCD coefficient {O^/'^^l^Si) = 1.12 GeV^ as the same as 
in Ref. [2^1 ■ Note, that this value is a little smaller (30 %) than the value which was used 
in our paper |22]- Concerning a charmed quark mass, the situation is not clear up to the 

M 

end. From one hand, in the nonrelativistic approximation one has = — , but there 

are many examples of taking smaller value of a c-quark mass in the amplitude of a hard 
process, for example, rric = 1.4 GeV. Taking into consideration above mentioned we perform 
calculations at rric = 1.5 GeV. The cinematic region under consideration is determined by 
the following conditions: < 1 GeV^, 60 < 1^ < 240 GeV, 0.3 < ^ < 0.9 and pr > 1 GeV, 
which correspond to the HI Collaboration data We assume that the contribution of the 
color octet mechanism is large at the z > 0.9 only. In the region of the small values of the 
z < 0.2 the contribution of the resolved photon processes [^0] as well as the charm excitation 
processes [21] may be large, too. All of these contributions are not in our consideration. 

Figs. 7-10 show our results which were obtained as in the conventional parton model as 
well as in the /c-r-factorization approach with the different parameterizations of the uninte- 
grated gluon distribution function. The dependence of the results on selection of a hard scale 
parameter /i is much less than the dependence on selection of a c-quark mass and selection of 
a parameterization. We put /i^ = + in a gluon distribution function and in a running 
constant as{fi'^). 

The count of a transverse momentum of the BFKL gluons in the /cT^-factorization ap- 
proach results in a flattening of the pr-spectrum of the J/ ip as contrasted by predictions of 
the parton model. For the first time this effect was indicated in the Ref. [22] , and later in 
the Ref. j2S]- Fig. 7 shows the result of our calculation for the pt spectrum of the J/ip 
mesons. Using the fcr-factorization approach we have obtained the harder p^-spectrum of 
the J/ip than has been predicted in the LO parton model. It is visible that at large values of 
Pt only the /c-r-factorization approach gives correct description of the data [211 • However, it 
is impossible to consider this visible effect as a direct indication on nontrivial developments 
of the small-x physics. In the article (2^1 was shown that the calculation in the NLO ap- 
proximation using colliniar parton model gives a harder pt spectrum of the J/ ip meson, too, 
which will agree with the data at the large pt- 

In the fcj^-factorization approach JB parameterization jH] gives p-r-spectrum, which is 
very close to experimental data. From the another hand in the case of JS parameterization 
jU] the additional i^'-factor approximately equal 2 is needed. 

The z spectra are shown in Fig. 8 at the various choice of the px cut: > 2, 4 and 6 
GeV, correspondingly. The relation between the theoretical predictions and experimental 
data is the same as in Fig. 7. The kr factorization approach give more correct description 
of the data especially at large value of z where the curve obtained in the collinear parton 
model tends to zero. 

Fig. 9 shows the dependence of the total J/ip photoproduction cross section on W at 
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0.3 < z < 0.8 and pt > 1 GeV. The shape of this dependence agrees well with the result ob- 
tained using JS parameterization 9 or KMR [lOj parameterization. However, the predicted 
absolute value of the cross section a^p is smaller by factor 2 than obtained data [21] • The 
results of calculation using JB |8J or GRV [7] parameterizations are larger and coincide with 
the data [221 better. 

As it was mentioned above, the main difference between the fcT-factorization approach 
and the conventional parton model is nontrivial polarization of the BFKL gluon. It is ob- 
vious, that such a spin condition of the initial gluon should result in observed spin effects 
during the birth of the polarized J/ ip meson. We have performed calculations for the spin pa- 
rameter a as a function z or pt in the conventional parton model and in the fc^-factorization 
approach: 

«(^) = = ^ (30) 

dz dz dpT dpT 

Here atot = ai + ax is the total J/ip production cross section, aL is the production cross 
section for the longitudinal polarized J/ip mesons, ax is the production cross section for 
the transverse polarized J/ip mesons. The parameter a controls the angle distribution for 
leptons in the decay J/ip —>■ in the J/ip meson rest frame: 

l + acos^{e). (31) 



dcos{6) 

The theoretical results for the parameter a{z) are very close to each other irrespective of 
the choice of an approach or a gluon parameterization |23j . 

For the parameter a{pT) we have found strongly opposite predictions in the parton model 
and in the fc^-factorization approach, as it is visible in Fig. 10. The parton model predicts 
that J/ip mesons should have transverse polarizations at the large pr (^(pt) = 0.6 at the 
Pt = 6 GeV), but /cr-factorization approach predicts that J/ip mesons should be longitudi- 
nally polarized ( q;(pt) = —0.4 at the pt = Q GeV). The experimental points lie in the range 
< Pt < 5 GeV and they have the large errors. However, it is visible that a^pr) decrease as 
Pt changes from 1 to 5 GeV. This fact coincide with theoretical prediction obtained in the 
/cT-factorization approach. Nowadays, a result of the NLO parton model calculation in the 
case of the polarized J/ip meson photoproduction is unknown. It should be an interesting 
subject of future investigations. If the count of the NLO corrections will not change predic- 
tions of the LO parton model for a{pT), the experimental measurement of this spin effect 
will be a direct signal about BFKL gluon dynamics. 

Nowadays, the experimental data on J/ ip polarization in photoproduction at large pt are 
absent. However there are similar data from CDF Collaboration [^21, where J/ip and ip' px- 
spectra and polarizations have been measured. Opposite the case of J/ip photoproduction, 
the hadroproduction data needs to take into account the large color-octet contribution in or- 
der to explain J /if) and ip' production at Tevatron in the conventional collinear parton model. 
The relative weight of color-octet contribution may be smaller if we use fc^-factorization ap- 
proach, as was shown recently in EH ISHl EH] • The predicted using collinear parton model 
transverse polarization of J/ip at large px is not supported by the CDF data, which can 
be roughly explained by the /c-r-factorization approach j^lj. In conclusion, the number of 
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theoretical uncertainties in the case of J/ip meson hadroproduction is much more than in 
the case of photoproduction and they need more comphcated investigation, which is why the 
future experimental analysis of J/ip photoproduction at THERA will be clean check of the 
coUinear parton model and the /cr-factorization approach. 
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Figure 1: The unintegrated gluon distribution function ^{x, k^, /j,'^) versus x at the fixed 
values of /i and k^. 
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Figure 2: The unintegrated gluon distribution function ^{x, k^, /i'^) versus k^ at the fixed 
values of and x. 
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Figure 3: The r) spectra of the D* meson at the various cut on a transverse momentum 
{Pt > 2, 4, 6 GeV, correspondingly from up to down) and 130 < W < 280 GeV. 
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Figure 5: The theoretical predictions for the W^-spectra at the various cut on the D* meson 
transverse momentum (p^ > 2,4, 6 GeV, correspondingly from up to down) and \r]\ < 1.5. 
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Figure 6: Diagrams used for description partonic process ■y + g —>■ J/ip + g. 
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Figure 7: The J/ijj spectrum on at the 60 < < 240 GeV and 0.3 < ^ < 0.9. 
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Figure 8: The J/^p spectrum on z at the 60 < < 240 GeV ( pr > 1, 2, 3 GeV, corres 
ingly from up to down). 
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Figure 9: The total J/%Ij photoproduction cross section versus W at the 0.3 < z < 0.8 and 
Pt > 1 GeV. 
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Figure 10: Parameter o; as a function of at the 0.3 < ^ < 0.9, 60 < < 240 GeV. 



18 



